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Signal Detection 

This invention relates to a method and apparatus for detecting signals. 
The invention is especially, but not exclusively, applicable to systems which 
5 may be subject to an unknown Doppler frequency shift affecting modulated 
coherent signals used for communication and ranging purposes, and 
particularly substantially continuous signals modulated by random or chaotic 
waveforms utilized by obstacle detection or collision avoidance systems. 

Such systems are designed to operate in multiuser and often hostile 

10 environments, and are intended for a wide range of applications, including 
automotive applications, industrial robotics and unmanned vehicle navigation. 

Fig. 1 is a block diagram of a typical microwave obstacle-detection 
system. The system comprises a signal generator 1 that produces a 
substantially continuous waveform x(t) with suitable bandwidth to provide 

15 required range resolution. The waveform x(t) may be deterministic (periodic 
or aperiodic), chaotic or purely random. 

The system also has a microwave oscillator 2 that generates a 
sinusoidal signal with required carrier frequency, a modulator 3 that 
modulates one or more of the parameters (such as amplitude, phase, or 

20 frequency) of the carrier signal with the modulating waveform x(t), 2l power 
amplifier (PA) 4 that amplifies the modulated carrier signal to a required 
level, a microwave transmit antenna (TA) 5 that radiates an electromagnetic 
wave representing the modulated carrier signal towards an obstacle 6, a 
microwave receive antenna (RA) 7 that receives an electromagnetic wave 

25 reflected back by the obstacle 6, an input amplifier (lA) 8 that amplifies a 
signal provided by the receive antenna (RA) 7, and a coherent demodulator 9 
that processes jointly the reference carrier signal supplied by the oscillator 2 
and a signal supplied by the input amplifier (lA) 8 to reconstruct a time- 
delayed replica y(t) of the modulating waveform x(t). 

30 The modulating waveform x(t) and its time-delayed replica y(t) are 

then processed jointly during a specified time interval by a suitable processor 
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10, such as correlator, to produce an estimate of the unknown time delay that 
is proportional to the distance (range) between the system and the obstacle 6. 

Fig. 2 shows an example of a correlation function of a synchronous 
random binary waveform. 
5 When there occurs a relative movement between a ranging system and 

an obstacle of interest, an electromagnetic wave reflected back by the obstacle 
and received by a coherent system will exhibit a Doppler frequency shift. The 
value codo of this (angular) frequency shift can be determined from: 

c 

10 where is the radial speed (i.e., the range rate) of a relative movement 

between the system and the obstacle and Oq is the (angular) carrier frequency 

of a transmitted electromagnetic wave having velocity c. 

The signal reflected back by a moving obstacle can be expressed as: 

Zr(0 = yxft - ro) cos [(CDo-^CODoXt- to) + ff] 
15 where y is the round-trip attenuation, tq is the delay corresponding to the range 
D, codo is the Doppler frequency shift, and ^ is an unknown constant phase 
shift. Strictly speaking, the value of tq cannot be constant for a nonzero 
Doppler frequency codo- However, in most practical cases it is assumed that 
the small changes in range D, of the order of the wavelength of the carrier 
20 frequency, cannot be discerned when determining a short-time time-delay 
estimate. Such assumption justifies decoupling delay and Doppler frequency 
measurements. 

A baseband signal corresponding to the received signal z,(t) has the 

form 

25 y(0 = x(t' To) cos(codo^-^(p) 

where ^ is an unknown phase shift. 

The correlator performs the following operation 

1 '? 

R:^ (^) = J X(.t - T)x(t - To ) COSiCOoot + <p)dt 

where the integral is evaluated for a plurality of hypothesized time delays 
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Xmin < T < Tmax- When the observation interval Tq is much shorter than one 
period (27c/fi)^o) of the Doppler frequency, the value of cos{Q}^Qt^+g>) is 
almost constant during Tq, so that the correlation integral can be 
approximated by 

1 

5 R^(r) - — cos(fi;^^f„ + (p)j x(t - r)x(t - r^)dt 

^0 0 

where the time instant is taken in the middle of the observation interval 7^ . 

When the correlation integral is calculated repeatedly for successive 
short processing intervals, each of duration 7J), the sequence of observed 
correlation functions may be represented by the plot of Fig. 3. The rate of 
10 change (in time) of the correlation function /?^(t) will correspond to the 

Doppler frequency co^^. The value of this frequency can be determined by 

applying some suitable form of spectral analysis. 

Fig. 4 is a block diagram of a conventional correlator, comprising 
variable delay line 1 1, multiplier 12 and integrator 13, followed by a spectrum 
15 analyser 14. When the total number of successive short processing intervals 
is large enough, the frequency spectram S(cd), observed at the output of the 
analyser, will exhibit a pronounced peak at the Doppler frequency (Oj^q . 

Fig. 5 is a block diagram of a multichannel correlator that uses a 
tapped delay line 15 to cover the entire interval of hypothesized time delays, 

20 rmin< T < tmax, in / stcps of Unit delay A using delay circuits 16. Multipliers 1 7 
process the delayed signals wiith y(t) and provide the outputs to integrators 
18. The required spectral analysis may be performed by a digital processor 19 
implementing the Discrete Fourier Transform (DFT). The principle of 
operation is thus similar to that of Fig. 4. 

25 The systems shown in Figs. 4 and 5, and also other similar known 

systems, attempt to approximate the correlation integral with two parameters 
of interest, delay time x and Doppler frequency cOj^ by combining a time 
sequence of correlation functions determined over a number of relatively 
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short observation intervals . Because of this approximation, such systems 

can only provide suboptimal solutions to the problem of joint estimation of 
time delay and Doppler frequency. 

Another prior art technique involves determining values of the 
5 following correlation integral 

1 

^xy ^z> ; 9^) = — J ^(^ - ^) cos(CD^t + q>)y{t)dt 
for the entire interval of hypothesized time delays 

min "^^ X < Tmax » 

and also for the entire interval of hypothesized Doppler frequencies 

10 ©Dmin < G>D < ©Dmax ; 

the unknown phase tp should also be varied over the interval (0, 2n). In this 
case, values of the correlation integral, determined for some prescribed range 
of argument values {%cod), define a two-dimensional correlation function. The 
specific values of arguments r and cooy say to and codoj that maximise the two- 
15 dimensional correlation function provide estimates of unknown time delay 
and unknown Doppler frequency. 

Fig. 6 is a block diagram of a suitable system capable of determining 
values of the correlation integral for some specified time delays r and Doppler 
frequencies cod* The system is formed by a correlator like that of Fig. 4, 
20 including variable delay line 11, multiplier 12 and integrator 13. However, 
the signal x(t) is first applied to a multiplier 20 which multiplies it by a 
variable Doppler signal cos(fi>^r + q>) . 

The procedure of adjusting the unknown value of phase q) can be 
avoided by suitably combining the values of the following two integrals: 

1 

25 /?^.,(r,fi>^;^) = — Jy(0^(r-T)cos(fi>^ryf (1) 

Rxys ; (P) = 7" I yit)x{t - r) sm{cDj,t)dt (2) 
evaluated for the entire interval of hypothesized time delays 
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and the entire interval of hypothesized Doppler frequencies 

0)Dmin < CDd < ©Dmax- 

In many industrial and automotive applications, systems capable of 
5 detecting moving obstacles and employing substantially continuous 
modulated microwave (or other coherent radiation), developed in accordance 
with the prior art, will be too complicated and also too expensive. 

Furthermore, obstacle-detection systems designed for multiuser 
environments will preferably exploit some suitable form of random or chaotic 
10 modulation which may preclude the use of conventional Doppler signal 
processors constracted in accordance with the prior art. 

Aspects of the present invention are set out in the accompanying 

claims. 

In accordance with a further aspect of the invention, representations of 
IS a transmitted signal and a received signal are processed together in such a 
manner as to determine whether or not a target signal is present (this target 
signal possibly representing an object reflecting the transmitted signal). An 
auxiliary representation, for example formed of one or more auxiliary signals, 
is introduced into the process. The auxiliary representation includes 
20 frequencies distributed throughout a range which corresponds to anticipated 
frequency modifications (for example Doppler shifts) in the target signal, as a 
result of which the system has response characteristics such that a target 
signal exhibiting frequency modifications within this range will give rise to a 
significant system output, whereas frequency modifications outside this range 
25 will not. 

Such an arrangement avoids the need for determining the value of the 
frequency modification (e.g. Doppler shift) but enables the rapid detection of 
any target signal incorporating a frequency modification of interest, using a 
simple stracture. The above-mentioned arrangement is thus particularly 
30 useful in applications where it is sufficient to discriminate between stationary 
obstacles and obstacles moving with specified velocities, and for applications 
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in which obstacles are placed into a number of classes according to obstacle 
velocity, for example zero/low, moderate/average, high/very high. 

Preferred aspects of the present invention are directed to the nature of 
the auxiliary representation and give rise to a number of advantages as 
5 explained below. 

Preferably, the auxiliary representation is in the form of finite-duration 
signal portions comprising multiple frequencies within the range of interest, 
and preferably distributed throughout this range. The dvu-ation of the signal 
portions is selected so as to maintain orthogonality of the components of the 

10 respective frequencies. By using finite-duration portions and selecting the 
shapes of these portions appropriately, it is possible to extend the system 
response characteristics such that a significant output is produced for not only 
target signal frequencies equal to those of the signal components used to 
constmct the auxiliary representation, but also adjacent frequencies. 

15 Preferably, each signal portion is created using a window function 

which results in a significant frequency response for all frequencies within the 
range of interest. However, it is difficult to arrange for a signal portion to 
produce a flat response throughout the frequency range of interest. 

For this reason, in accordance with a preferred embodiment of the 

20 invention, the auxiliary representation comprises different finite-duration 
signal portions, each comprising components having a respective set of 
frequencies, and arranged such that they give rise to complementary 
frequency response characteristics. Accordingly, by combining the results 
achieved using the respective different auxiliary signal portions, it is possible 

25 to get an overall frequency response characteristic which is significantly 
flatter than if only one type of signal portion were used. 

Preferably, each signal portion comprises components of frequencies 
interleaved with those of the other signal portion or portions. 

An indicated by Fig. 3, at any given instant the correlation function 

30 may have a low level depending upon the phase of the Doppler shift. In order 
to ensure a significant output irrespective of the unknown phase, the auxiliary 
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representation preferably comprises a first representation and a second 
representation, the components of respective frequencies in the first 
representation having a quadrature phase relationship with the corresponding 
frequency components in the second representation. By combining the 
5 responses achieved as a result of theses two auxiliary signals, the output is 
substantially unaffected by the phase of the Doppler shift. 

Accordingly, in the preferred embodiment, there are at least four types 
of auxiliary signal portions: signal portions which differ in that they comprise 
frequencies interleaved with frequencies of other portions, and signal portions 

10 which differ in that their frequency components are in quadrature relationship 
with those of other portions. 

The preferred embodiment uses four types of signal portions, but it 
would be possible to produce a system which uses more types, while 
maintaining the advantages indicated above. For example, the frequency 

15 components may be split into three or more sub-sets. 

It would be desirable to synthesise digitally the auxiliary 
representation. In order to facilitate this, the auxiliary representation 
preferably exhibits small values of the peak factor, i.e., the peak value divided 
by its root-mean square value. Minimising this enables the waveform to be 

20 represented digitally with a minimum number of bits for a specified 
reconstruction error. 

A small peak factor can be achieved by appropriate selection of the 
initial phase values of the signal components for each signal portion. 

In the preferred embodiment described below, the transmitted and 

25 received representations are signals which are combined by correlation. 
However, other arrangements are possible. For example, the received signal 
may be modified (e.g. multiplied) by the auxiliary representation and then 
delivered to a filter matched to the characteristics of the transmitted signal. 
This use of matched filters, which utilise a representation of the transmitted 

30 signal rather than the signal itself, is, per se, known in the art. 
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In another alternative embodiment, the correlator is replaced by a 
suitably modified version of the time-delay discriminator described in WO-A- 
(X)/39643, the contents of which are incorporated herein. 

Although the present invention is being described primarily in the 
5 context of detection of obstacles which produce a Doppler shift in a reflected 
signal, the invention is applicable also in other areas. For example, the 
invention is applicable to communication systems in which the local 
oscillators in the transmitter and receiver may not be accurately matched; a 
frequency offset between the oscillators would give rise to a frequency shift 
10 similar to a Doppler shift. The invention thus permits use of equipment which 
has wider tolerances, and is thus less expensive, than prior art arrangements. 

The invention relates both to a method of signal processing and to 
apparatus operating according to such a method. 

Arrangements embodying the invention will now be described by way 
15 of example with reference to the accompanying drawings. 

Fig. I is a block diagram of a typical microwave obstacle-detection 

system. 

Fig. 2 is a correlation function of a synchronous random binary 
waveform. 

20 Fig. 3 shows the effects of different values of phase of the Doppler 

shift in a target signal on the correlation function of a synchronous random 
binary waveform. 

Fig. 4 is a block diagram of a conventional correlator followed by a 
spectrum analyser. 

25 Fig. 5 A block diagram of a multichannel correlator followed by a 

digital processor implementing the Discrete Fourier Transform. 

Fig. 6 is a block diagram of a system for determining values of the 
correlation integral for various time delays and various Doppler frequencies. 

Fig. 7 is a representation of the time-frequency plane divided into 
30 range-Doppler resolution cells. 
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Fig. 8 is a representation of the time-frequency plane divided into 
range resolution cells for a specified interval of Doppler frequencies. 

Fig. 9 is a block diagram of a two-channel Doppler-tolerant signal 
processor in accordance with the present invention. 

Fig. 10 is a block diagram of a multichannel correlator constructed in 
accordance with the prior art. 

Fig. 11 is an example of an auxiliary representation in the form of 
signal portions for use in the processor of Fig. 9. 

Fig. 12 schematically illustrates an example of operations performed 
by a channel combiner of the processor of Fig. 9. 

Fig. 13 shows an alternative configuration of a two-channel Doppler- 
tolerant signal processor in accordance with the present invention. 

Fig. 14 is a block diagram of a four-channel Doppler-tolerant signal 
processor in accordance with the present invention. 

Fig. 15 is a block diagram of a quadrature Doppler generator for use in 
a Doppler-tolerant signal processor in accordance with the present invention- 
Figs. 16 (a) and (b) and Figs. 17 (a) and (b) show examples of 
synthesized signal portions used to construct the auxiliary representation. 

Fig. 18 shows the partial frequency response of a two^channel 
Doppler-tolerant signal processor employing the signal portions of Figs. 16 
(a) and (b). 

Fig. 19 shows the complementary frequency response of a two- 
channel Doppler-tolerant signal processor employing the signal portions of 
Figs. 17 (a) and (b). 

Fig. 20 shows the combined frequency response of a two-channel 
Doppler-tolerant signal processor employing all four signal portions. 

Fig. 21 shows ripples of the combined frequency response of Fig. 20. 

Fig. 22 shows the combined frequency response of a two-channel 
Doppler-tolerant signal processor employing a squared cosine window. 

Fig. 23 shows ripples of the combined frequency response of Fig. 22. 
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Fig. 24 shows a modified Kaiser window (solid line) and a squared 
cosine window (broken line). 

Fig. 25 shows two sets, A and B, of auxiliary signals employed by two 
two-channel Doppler-tolerant signal processors operating in a tandem 
configuration. 

Fig. 7 shows a representation of the time-frequency plane divided into 
range-Doppler resolution cells and a hypothetical obstacle occupying a single 
cell. Prior art arrangements are designed to establish the coordinates to,codo of 
the cell. 

When obstacles are to be classified into only two alternative classes: 
obstacles movingj vith spe cified velocities and all other obstacles, the time- 
frequency plane wiU-b^-divided^i^to a numb erof re solution cells equal to the 
number of requi red range cells. An obstacle det ectedTlln^^y of those 
resolution cells will have the range specifie d bv that celL but its velocity may 
be one of the plurality of admissible velocities^QQcresponding to an entire 
interval of Dopplei;;i&equ6B64es»^^ay, coomin < coo < coomax^ 

Fig. 8 is a representation of the time-frequency plane divided into a 
number of range resolution cells for a specified interval of Doppler 
frequencies. An obstacle located at the range cell to may have a radial velocity 
within an interval of velocities corresponding to the entire Doppler interval 

CODmin <COd< CODmax- 

In order to detect obstacles moving with any of the specified 
velocities, the two correlation integrals, given by (1) and (2), can be replaced 
by 

1 

1 ^° 

Rrys ^2> ; 9>) = — J yit)x{t - r)ds {t)dt (4) 

•'o 0 

where dc(t) and ds(t) are especially designed auxiliary signals 
comprising a plurality of suitable frequency components. 
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The auxiliary signals dc(t) and ds(t) should be designed in such a way 

as: 

- to contain all, mutually orthogonal within the interval To, frequency 
components covering the entire interval of Doppler frequencies of interest; 

- to contain two quadrature components at each frequency of interest, 
thus providing independence of unknown phase shift; 

- to exhibit small values of the peak factor (i.e., the peak value of a 
waveform divided by its root-mean-square value) required to represent an 
analogue waveform digitally with a minimum number of bits for a specified 
value of a reconstraction error; 

- to provide a signal processor with a flat frequency response over the 
entire interval of frequencies of interest, and also a required attenuation at all 
other Doppler firequencies. 

According to a preferred embodiment of the present invention, the two 
auxiliary signals, dc(t) and ds(t), are obtained by suitably combining four 
finite-duration chirps, dcE^O , dsE(t) , dco(t) and dsd(t), defined by: 

K 

dcE (0 = {t)Y,oos[27ri2k + «)/^/ + ] (5) 

K 

dsE (0 = -W£ (02 sin[2;r(2fc + w)/pr + aj (6) 
dco (0 = (0i;cos[2;r(2fc + 1 + + A 1 (7) 

dsoiO=-Wo(t)f,sm[27Ci2k + l+u)f^t + ^,] (8) 

where w^it) and w^(t)arc suitable window (or tapering) functions, K is the 
number of even or odd frequency components, u is the integer frequency shift, 
is the fundamental frequency, and cTj^ and yffj^ are phase angles so 

selected as to obtain the minimum, or at least low, value of the peak factor. 

The chirp dsE(t) is a quadrature copy of the corresponding chirp dcdt)\ 
similarly, the chirp dso(t) is a quadrature copy of the corresponding chirp 
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dco(t)' All harmonics of the fundamental frequency are shifted (by u) and 

split into two sets, each comprising either even or odd harmonics. Those sets 
are then employed to synthesize two orthogonal groups of chirps: {rfc£(0» 
dsE{t)}. and {dco(t). dso(t)}. 

Window functions w^(t)and w^C/) belong to a large class of known 

and widely used window functions, such as Blackman windows, Kaiser- 
Bessel windows or Dolph-Chebyshev windows. The duration Tw of a 
selected window is so chosen as to exploit the orthogonality of the frequency 
components, and the window type and shape are mainly determined by a 
required frequency resolution at each Doppler frequency of interest. 

In order to simplify the optimisation procedure and also to obtain 
symmetrical chirps, the values of phases oTj^ and J3f^ can be restricted to only 

two values, 0 and tu. In such a case, the four chirps, dcE^t) , dsE(t) , dco(t) and 
ds(y(t) will be synthesized as follows: 

K 

dcEiO = w^COX^* cos[2;r(2fc + u)//] (9) 

dsE (0 = -^E (OX ^* sin[2;r(2A: + m)/^?] (10) 

dcoit) = Wo it)i,K oos{2jt{2k + 1 + (1 1) 

dso it) = -Wo (0£ sin[2;r(2fc + 1 + «)/^/] (12) 

where the coefficients and b^^ may assume only two values +1 or 

-1, and are so selected as to minimize the value of the peak factor. 

A signal processor determining and combining the values of the two 
correlation integrals (3) and (4) will be referred to as a "Doppler-tolerant" 
signal processor. The resulting wideband Doppler compensation will make 
all obstacles of interest detectable ( Visible 0 at the output of a Doppler-tolerant 
signal processor, while all other obstacles will be strongly attenuated. 
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Fig. 9 is a block diagram of a two-channel Doppler-tolerant signal 
processor, incorporating a quadrature Doppler generator QDG, according to 
the present invention. The processor utilizes two identical signal processing 
blocks, CRI and CRQ, either of which may be a conventional correlator or 
5 another appropriate signal processor. In the following, it will be assumed that 
each of the blocks CRI and CRQ is a conventional multichannel correlator. 

Fig. 10 is a block diagram of a multichannel correlator, constructed in 
accordance with the prior art (see Fig. 5), which is an example of a suitable 
structure for each of blocks CRI and CRQ. The correlator provides J discrete 
10 values of the correlation function detemiined at the multiples of unit delay A. 

In the two-channel Doppler-tolerant signal processor of Fig. 9, the 
function of Doppler compensation is implemented jointly by a quadrature 
Doppler generator QDG, two analogue multipliers, AMI and AMQ, and a 
channel combiner IQC. 
15 A complete set of operations is performed by the two-channel signal 

processor during a 'basic processing time interval 27V' split into two equal 
subintervals T^, referred to as an even cycle and an odd cycle. 

A baseband signal y{t), supplied by the coherent demodulator 9 in Fig. 
1 (shown here as comprising a mixer AM for mixing the input signal with a 
20 local oscillator signal, followed by a low pass filter LPF), is applied 
simultaneously to inputs II and Ql of the multipliers AMI and AMQ, 
respectively. The other two inputs, 12 and Q2, of the multipliers are driven by 
auxiliary signals produced by the quadrature Doppler generator QDG. 

During an even cycle, the generator QDG produces simultaneously 
25 two chirps: dc^t) at output 12, and ds^t) at output Q2. During an odd cycle, 
another two chirps, dccAt) and dsoUh are produced simultaneously at outputs 
12 and Q2, respectively. When the processing interval 2Tw is replicated, the 
alternate cycles contain chirps dcd^) and dco(t) which will thus form a 
periodic auxiliary signal dc(t). Similarly, a periodic auxiliary signal d^t) will 
30 be obtained by periodic extension of combined chirps dsEi(t) and dso(t). Each 
of these two resulting quadrature signals, dc(t) and ds(t)^ contains the same 
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frequency components covering the entire hypothesized Doppler frequency 
range of interest. 

The generator QDG also produces a pulse train EOC to indicate the 
time instants of cycle transitions, but not necessarily a particular even or odd 
5 cycle. 

Fig. 11 shows an example of two auxiliary quadrature signals, dc(t) 
and ds(t), and a pulse train EOC produced by the quadrature Doppler 
generator QDG. 

Each correlator, CRI or CRQ, processes jointly a reference signal x(t) 
10 and a signal supplied by the corresponding multiplier, AMI or AMQ, during 
even and odd cycles. These cycles, each of duration TV, are initiated by 
consecutive pulses of the pulse train EOC. Results presented at the output of 
the correlator CRI, after completing an even cycle and an odd cycle, are 
denoted by CE and CO, respectively. Similarly, SE and SO denote results 
15 obtained at the output of the correlator CRQ after completing an even cycle 
and an odd cycle, respectively 

During an even cycle, the channel combiner IQC determines a 
representation of 

RE = ylCE''-hSE^ 
20 and during an odd cycle, IQC determines a representation of 

The above representations can be determined with the use of a 
conventional Pythagoras processor' or another suitable processor constmcted 
according to the prior art. 
25 In order to simplify the structure of the channel combiner IQC, the 

values of RE and RO can also be determined from some suitable 
approximations. For example, 

RE ^ max{|CE|,|SE|} + 0.5 mm{|CE|,|SE|} 
RO - max{|CO|,|SO|} + 0.5 min{|COUSO|} 
30 Other approximations can also be constructed as known in the art. 
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Representations of RE and RO, determined for a prescribed number L 
of the processing intervals 27^, are suitably combined together, e.g., added or 
averaged, and then the resulting combined representation REO is provided at 
the output of the combiner IQC. For example, for a single observation 
interval, L = 1 and REO = RE + RO. All operations performed by the channel 
combiner can be implemented either by a special-purpose digital processor or 
by a standard programmable digital signal processor. 

For obstacle detection purposes, the combined representation REO is 
to be compared to a predetermined decision threshold. Obstacle detection 
will be declared in those delay (range) cells, where the threshold has been 
exceeded by the corresponding value of the combined representation REO. 

Fig. 12 is a schematic illustration of an example of operations 
performed by the channel combiner IQC and the resulting detection decisions, 
which employ a decision threshold DT. 

Fig. 13 shows an alternative configuration of a two-channel Doppler- 
tolerant signal processor with Doppler compensation. In this case, it is the 
reference signal x(t)^ and not a demodulated signal y(t)^ that is multiplied by 
the chirps: dcdt)^ ds^th dco(0 and dsoO) produced by the quadrature Doppler 
generator QDG. Except for this modification, all operations and functions 
performed by the processor remain identical to those performed by the 
configuration shown in Fig. 9. 

Fig. 14 is a block diagram of yet another configuration of a Doppler- 
tolerant signal processor utilizing Doppler compensation. The system 
comprises four channels: a pair of in-phase channels and a pair of quadrature 
channels, each with their own correlator CRBE, CRQE, CRIO or CRQO. In 
this case, the chirps: dcE{t), dsdth dco(t) and dso(t) are supplied 
simultaneously during a single cycle T^. Consequently, the quadrature 
Doppler generator has now four signal outputs, in addition to an output for a 
processing interval signal PI indicative of the interval over which correlation 
should take place. Using additional separate channels for even-cycle chirps 
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and odd-cycle chirps reduces by half the time required for determining the 
combined representation REO by a suitably modified channel combiner IQC. 

It is also possible to devise another configuration of the quad-channel 
Doppler-tolerant signal processor processor shown in Fig. 14, by exploiting 
5 the modification presented in Fig. 13. 

In some applications, it may be convenient to employ a single-channel 
Doppler tolerant signal processor, in which the four chirps, dcE(t)y dsd^)y 
dco(0 and dso(t) will be produced and used sequentially. 

All the above configurations of a Doppler-tolerant signal processor, 

10 and also other similar structures, will utilize an appropriate form of a 
quadrature Doppler generator QDG to provide suitable auxiliary signals for 
Doppler frequency compensation in accordance with the present invention. 

Fig. IS is a block diagram of a suitable quadrature Doppler generator 
QDG. The generator comprises an JV-bit binary counter BCT, driven by a 

15 clock generator CLK, two memories, CVM and SVM, two buffer registers, 
CBR and SBR, two digital-to-analogue converters, CDA and SDA, and a 
trigger circuit TRG. ^ 

The operations performed by the quadrature Doppler generator can be 
summarised as follows. The most significant bit (MSB) of the output of the 

20 free-running counter BCT identifies even cycles, when MSB = 0, and odd 
cycles, when MSB = 1. The MSB and remaining (AT-l) bits are used jointly to 
define addresses of 2^ cells of each memory, CVM and SVM. The memory 
CVM contains digital amplitude samples of a chirp dcE(0 at addresses with 
MSB = 0, whereas samples of a chirp dco(0 are stored at addresses with MSB 

25 =1. Similarly, the memory SVM contains digital samples of a chirp dsE(0 at 
addresses with MSB = 0, whereas samples of a chirp dso(0 are stored at 
addresses with MSB = 1. Consequently, each of the four chirps is represented 
by 2^'^ binary words, each comprising M bits, where the number Af is chosen 
to attain a specified accuracy of the resulting digital representation of chirp 

30 samples. 
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Successive digital words appearing at the outputs of memories CVM 
and S VM are transferred to the corresponding buffer registers, CBR and SBR, 
at the time instants determined by a suitably delayed pulse train available at 
output 2 of the clock generator CLK. Digital words held in the buffer 
5 registers CBR and SBR are then converted respectively into analogue samples 
by two corresponding digital-to-analogue converters, CDA and SDA. As a 
result, the converter CDA produces repeatedly: 2^"^ analogue samples of chirp 
dcdt) followed by 2^"' analogue samples of chirp dco(t)* Similarly, the 
converter SDA produces repeatedly: 2^"^ analogue samples of chirp dsE/lO 

10 followed by 2^'^ analogue samples of chirp dso(0 • 

A trigger circuit, TRG, generates a short pulse at the time instants 
coinciding with each transition of a binary waveform representing values of 
the most significant bit (MSB) of successive states of the counter BCT. A 
sequence of such short pulses, denoted by EOC, is used to indicate the 

15 beginning of each (even or odd) cycle. The pulse train EOC is utilized in the 
Doppler-tolerant signal processor, shown in Fig. 9 and Fig. 13, to define 
integration intervals of the correlators, CRI and CRQ, and also to initiate 
various operations performed by the channel combiner IQC. 

For a quad-channel Doppler-tolerant signal processor, the quadrature 

20 Doppler generator will be modified as follows. The number of bits of the 
binary counter BCT will be reduced to (A^-1), and successive counter states 
will determine addresses of four memories containing digital samples of 
corresponding four chirps, dcE(0 , dsE(0 , dco(0 and dso(t). The output of 
each of the four memories will then drive a separate digital-to-analogue 

25 converter via a suitable buffer register. As a result, analogue representations 
of the four chirps will be available simultaneously at four outputs of such 
modified quadrature Doppler generator. 

The modifications of the quadrature Doppler generator required for a 
single-channel operation will be obvious to those skilled in the art. 

30 In order to illustrate some basic aspects of the present invention, and 

also to demonstrate its potential benefits, consider the following example. 
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Let the number of frequency components jfir=31 and the integer 
frequency shift w=4; hence, the minimum frequency and the maximum 
frequency are equal to 6fp and 61 fp, respectively. Assume that each of the 
four chirps is to be represented by 2048 discrete samples and that w^t) = 
wo(t) = Mty> assume also that a^^ = b^^ for fe=l,2,. . ..,31. 

It has been found that the peak factor achieves its minimum value, 

when 

ai=+l ^2 = +! a3 = +l a4=+l 05 = +1 a6 = +l ci? = +1 
fl8 = -l a9 = -l aio = -l aii=+l ai2 = +l ^13 = +1 
au = -1 ^15 = -1 ^16 = -1 «17 = +1 ^18 = +1 ^19 = -1 

= +1 CL2\ = +1 ^22 = -1 ^23= "1 ^24 = +1 ^25 = "1 
026 = -1 a27 = +l 028=-! ^29=+! 030 = -! a31=+l- 

Application of the above coefficients has resulted in the reduction of 
the peak factor value from 31/V(31*0.5) ^ 7.9 (when all components are 
simply added with unit coefficients) to only 1.9. It is desirable to ensure the 
peak factor value is less than 2. 

To achieve a flat frequency response over the selected frequency 
range, and provide a required attenuation at all other frequencies, the discrete 
window function w(j) is a modified Kaiser window of the form 

/o(^)-l 

where m = (j - 1024)/1024 for y=0,l ,2047, /o(») is the modified Bessel 

function, and the shape parameter ^ is equal to 6.24. 

Samples of the four discrete chirps, dcElU) . dssU) . dcoU) dsoO). 
are determined from 

dcEU) = 2wU)^a, cos[2;r(2A: + 4)j71024] 
dsE U) = -2w(j)X«t sin[2;r(2A: + 4)j / 1024] 
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dcoU) = 2wU)^a, cos[2;r(2A: + 5)y71024] 

dsoU) = -2wU)^a, sin[2;r(2/: + 5);71024] 

for7=0,l,....,2047. 

The chirps rfcfO") » ^^fOV » ^coC/V and dsdU), synthesized above, can be 
produced by a quadrature Doppler generator QDG, with parameters ^=12 and 
M=12. 

The synthesized chirps dcEij) , dssij) , dcoQ) and dsoU) are shown in 
Fig. 16 and Fig. 17. As seen, the chirps exhibit various forms of symmetry 
which can be exploited to reduce the memory of the quadrature Doppler 
generator QDG, at the expense of more complicated timing and control 
functions. 

The chirps dcsij) and dsE/U) are used for correlation during an even 
cycle to produce the results CE and SE, respectively. Those results are then 
combined by the channel combiner IQC as follows 

The resulting partial frequency response of a two-channel Doppler- 
tolerant signal processor is shown in Fig. 18. 

Similarly, the chirps dcoij) and dsoU) are used for correlation during 
an odd cycle to produce the results CO and SO, respectively. Those results are 
combined by the chaimel combiner IQC as follows 

The resulting complementary frequency response of a two-channel 
Doppler-tolerant signal processor is shown in Fig. 19. 

The two frequency response characteristics are superimposed (added) 
to produce the combined frequency response, shown in Fig. 20. As seen, the 
relative attenuation at out-of-band frequencies can exceed 50 dB. 
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The plot of the combined frequency response in Fig. 21 reveals that 
the ripples of the response do not exceed ±0.25 dB at all frequencies of 
interest. 

Fig. 22 shows the combined frequency response when the modified 
5 Kaiser window is replaced a simple cosine window of the form 

2r j-1024 ;r^ 

>V( 7) = COS 

1^ 1024 2 J 

for y=0,l,.. ..,2047. 

The cosine window provides greater attenuation of stationary 

obstacles (when Doppler frequency is zero); however, the frequency response 
10 is not as flat as that corresponding to the modified Kaiser window. The shape 

and location of ripples are shown in Fig. 23. 

In all the above configurations of a Doppler-tolerant signal processor, 

the received signal y{t) is multiplied by chirps with nonimiform envelopes 

determined by a selected window. Because of this loss of signal energy, the 
15 detection performance of the system can be significantly reduced. The 

examination of various windows utilized in a Doppler-tolerant signal 

processor shows that their shape is close to that of a squared cosine function. 

For example. Fig. 24 shows that a modified Kaiser window and a squared 

cosine window almost overlap. This observation and the trigonometric 
20 identity sxx^tp + cos^^= 1 suggest that the loss due to the use of windows can 

be almost eliminated when a received signal is processed simultaneously by 

two Doppler-tolerant signal processors operating in tandem configuration. 

While either processor employs the same set of chirps, all chirps of a 

duplicated set are shifted in time by TJ2^ where is the window duration. 
25 Fig. 25 shows two sets of chirps employed by a tandem arrangement 

comprising two two-channel Doppler-tolerant signal processors, each 

arranged in accordance with Fig. 9. 

A more efficient system will utilize a tandem configuration 

comprising two quad-channel Doppler-tolerant signal processors, each 
30 developed in accordance with Fig. 14. 
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This technique, whereby auxiliary signals are used to investigate a 
primary signal, and wherein each auxiliary signal comprises finite-duration 
portions overlapping portions of the other auxiliary signal, is considered 
independently inventive. It can be applied to other areas of signal processing 
5 beyond the specific areas disclosed herein, including but not limited to 
arrangements in which the auxiliary signals contain components of different 
frequencies. It can also be extended to cover the use of more than two trains 
of auxiliary signals. ^ 

The correlators of the embodiments described above could be replaced 
10 by other arrangements such as a suitably modified time-delay discriminator as 
described in WO-A-00/39643. 



